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Abstract 



The first measurements of the invariant differential cross sections of inclusive it and 77 meson pro- 
duction at mid-rapidity in proton-proton collisions at y/s = 0.9 TeV and y/s = 7 TeV are reported. 
The it measurement covers the ranges 0.4 < p t < 7 GeV/c and 0.3 < p t < 25 GeV/c for these 
two energies, respectively. The production of 7} mesons was measured at y 7 ! = 7 TeV in the range 
0.4 < p t < 15 GeV/c. Next-to-Leading Order perturbative QCD calculations, which are consistent 
with the n° spectrum at y/s = 0.9 TeV, overestimate those of ft and rj mesons at y/s = 7 TeV, but 
agree with the measured rj /n° ratio at yfs = 7 TeV. 
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1 Introduction 

Hadron production measurements in proton-proton collisions at the Large Hadron Collider (LHC) Q] 
energies open a new, previously unexplored domain in particle physics, which allows validation of the 
predictive power of Quantum Chromo Dynamics (QCD) [2]. A quantitative description of hard pro- 
cesses is provided by perturbative QCD (pQCD) supplemented with parton distribution functions (PDF) 
f(x) and fragmentation functions (FF) D(z), where x is the fraction of the proton longitudinal momen- 
tum carried by a parton and z is the ratio of the observed hadron momentum to the final-state parton 
momentum. Due to the higher collision energy at the LHC, the PDF and FF can be probed at lower 
values of x and z, respectively, than in previous experiments. Such measurements can provide further 
constraints on these functions, which are crucial for pQCD predictions for LHC energies. Furthermore, 
while pion production at the Relativistic Heavy Ion Collider (RHIC) Q is considered to be dominated 
by gluon fragmentation only for p t < 5 — 8 GeV/c H |5j, at LHC energies it should remain dominant 
for p t < 100 GeV/c (HIT]]. Theoretical estimates [6] suggest that the fraction of pions originating from 
gluon fragmentation remains above 75 % in the p t range up to 30 GeV/c. Here, the measurement of the 
71° production cross section at LHC energies provides constraints on the gluon to pion fragmentation |[H 
in a new energy regime. In addition, the strange quark content of the v\ meson makes the comparison to 
pQCD relevant for possible differences of fragmentation functions with and without strange quarks l9l . 
Furthermore, the precise measurement of 71° and tj meson spectra over a large p t range is a prerequisite 
for understanding the decay photon (electron) background for a direct photon (charm and beauty) mea- 
surement. Finally, a significant fraction of hadrons at low p t is produced in pp collisions via soft parton 
interactions, which cannot be well described within the framework of pQCD. In this kinematic region 
commonly used event generators like PYTHIA iflOl or PHOJET ifTTIl have to resort to phenomenological 
models tuned to available experimental data delivered by lower-energy colliders like SppS, RHIC, and 
Tevatron |[T2l . to adequately describe hadron production. The large increase in center-of-mass energy at 
the LHC provides the possibility for a stringent test of the extrapolations based on these models. 

This paper presents the first measurement of neutral pion and v\ meson production in proton-proton 
collisions at center-of-mass energies of y/s = 0.9 TeV and 7 TeV in a wide p t range with the ALICE 
detector lPT3l . The paper is organized as follows: description of the subdetectors used for these mea- 
surements, followed by the details about the data sample, as well as about event selection and photon 
identification, is given in section[2j Section[3]describes the algorithms of neutral meson extraction, meth- 
ods of production spectra measurement, and shows the systematic uncertainty estimation. Results and 
their comparison with pQCD calculations are given in section [4] 

2 Detector description and event selection 

Neutral pions and r\ mesons are measured in ALICE via the two-photon decay channel. The photons 
are detected with two methods in two independent subsystems, with the Photon Spectrometer (PHOS) 
Ifl4l and with the photon conversion method (PCM) in the central tracking system employing the Inner 
Tracking System (ITS) [15] and the Time Projection Chamber (TPC) [16]. The latter reconstructs and 
identifies photons converted to e + e~ pairs in the material of the inner detectors. The simultaneous 
measurements with both methods with completely different systematic uncertainties and with momentum 
resolutions having opposite dependence on momentum provide a consistency check of the final result. 

The PHOS detector consists at present of three modules installed at a distance of 4.60 m from the in- 
teraction point. PHOS covers the acceptance of 260° < (p < 320° in azimuthal angle and |t]| < 0.13 in 
pseudorapidity. Each module has 3584 detection channels in a matrix of 64 x 56 cells. Each detection 
channel consists of a lead tungstate, PbWO/t, crystal of 2.2 x 2.2 cm 2 cross section and 18 cm length, 
coupled to an avalanche photo diode and a low-noise charge- sensitive preamplifier. PHOS operates at 
a temperature of —25 °C at which the light yield of the PbW04 crystal is increased by about a factor 3 
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compared to room temperature. PHOS was calibrated in-situ by equalizing mean deposited energies in 
each channel using events with pp collisions. 

The Inner Tracking System (ITS) iPTBl consists of six layers equipped with Silicon Pixel Detectors (SPD) 
positioned at a radial distance of 3.9 cm and 7.6 cm, Silicon Drift Detectors (SDD) at 15.0 cm and 
23.9 cm, and Silicon Strip Detectors (SSD) at 38.0 cm and 43.0 cm. The two innermost layers cover a 
pseudorapidity range of |r]| < 2 and |tj| < 1.4, respectively. 

The Time Projection Chamber (TPC) lfT6l is a large (85 m 3 ) cylindrical drift detector filled with a 
Ne/C02/N2 (85.7/9.5/4.8%) gas mixture. It is the main tracking system of the Central Barrel. For 
the maximum track length of 159 clusters it covers a pseudorapidity range of |tj| < 0.9 over the full 
azimuthal angle. In addition, it provides particle identification via the measurement of the specific ioni- 
sation energy loss (dE/dx) with a resolution of 5.5% ifTol . The ITS and the TPC are aligned with respect 
to each other to the level of few hundred using cosmic-ray and proton-proton collision data lfl5l . 

The event selection was performed with the VZERO detector [ 17"] in addition to the SPD. The VZERO 
is a forward scintillator hodoscope with two segmented counters located at 3.3 m and —0.9 m from the 
interaction point. They cover the pseudorapidity ranges 2.8 < r\ < 5.1 and —3.7 < tj < —1.7, respec- 
tively. 

The proton-proton collision data used in this analysis were collected by the ALICE experiment in 2010 
with the minimum bias trigger MBor lTT8l . This trigger required the crossing of two filled bunches and a 
signal in at least one of the two SPD pixel layers or in one of the VZERO counters. An offline selection 
based on time and amplitude signals of the VZERO detectors and the SPD was applied to reject beam- 
induced and noise background iTTSl . Pileup collision events were identified imposing a criterion based on 
multiple primary vertices reconstructed with the SPD detector, and removed from the further analysis. 
The cross sections for the MBor trigger have been calculated from other measured cross sections at the 
same energies with appropriate scaling factors. At y/s = 7 TeV the cross section for the coincidence 
between signals in the two VZERO detectors, <7mb A nd' was measured in a Van-der-Meer scan [19], and 
the relative factor cTmb and /<7mb r = 0.873 with negligible error as obtained from data was used. At 
■*/s = 0.9 TeV the cross section <Tmb r nas been calculated from the inelastic cross section measured in 
pp collisions at y/s = 0.9 TeV [20] and relative factor cJmb r /oinei = 0.91+qqi estimated from Monte 
Carlo simulations [l~9l . 

Table[T]shows the values of the cross section obtained at both energies as well as the integrated luminosity 
of the total data samples used. In the photon conversion analysis, only events with a reconstructed vertex 
(~90% of the total) are inspected, and those events with a longitudinal distance (i.e. along the beam 
direction) between the position of the primary vertex and the geometrical center of the apparatus larger 
than 10 cm are discarded. The analysis using PHOS as well as Monte Carlo simulations show that the 
number of % s in events without a reconstructed vertex is below 1% of the total number of 7T°s. 



V^(TeV) 


<*mBor ( mb ) 


a^ L (mb) 


0.9 


47.8+^(sysf) 


52.5±2(syst) 


7 


62.2±2.2(syst) 


73.2 ^ ±2.6 lumi 


\/s(TeV) 


££ (nb- 1 ) 


PCM 


PHOS 7T° 


PHOS tj 


0.9 


0.14 


0.14 




7 


5.6 


4.0 


5.7 



Table 1: Cross sections of the reactions and integrated luminosities of the measured data samples for the two beam 
energies (top), and luminosities used in the different analyses for the 7 TeV data (bottom). 
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To maximize the pion reconstruction efficiency in PHOS, only relatively loose cuts on the clusters (group 
of crystals with deposited energy and common edges) were used: the cluster energy was required to be 
above the minimum ionizing energy Cluster > 0.3 GeV and the minimum number of crystals in a cluster 
was three to reduce the contribution of non-photon clusters. 

Candidate track pairs for photon conversions were reconstructed using a secondary vertex (V0) finding 
algorithm [21 J. In order to select photons among all secondary vertices (mainly y, K^, A and A), electron 
selection and pion rejection cuts were applied. The main particle identification (PID) selection used the 
specific energy loss in the TPC (dE/dx). The measured dE/dx of electrons was required to lie in the 
interval [— 4a d£ / iY , +5a d£ /^J around the expected value. In addition, pion contamination was further re- 
duced by a cut of 2a above the nominal pion dE/dx in the momentum range of 0.25 GeV/c to 3.5 GeV/c 
and a cut of 0.5a at higher momenta. For the 7 reconstruction constraints on the reconstructed photon 
mass and on the opening angle between the reconstructed photon momentum vector and the vector join- 
ing the collision vertex and the conversion point were applied. These constraints were implemented as 
a cut on the X 2 (y) defined using a reconstruction package for fitting decay particles (22). The photon 
measurement contains information on the direction which allows to reduce the contamination from sec- 
ondary neutral pions. With the conversion method a precise y-ray tomograph of the ALICE experiment 
has been obtained [23]. The integrated material budget for r < 180 cm and |t]| <0.9 is 11.4±0.5% Xo as 
extracted from detailed comparisons between the measured thickness and its implementation in Monte 
Carlo simulations based on the GEANT 3.21 package using the same simulation runs for the material 
studies as for the 71° measurement. Photon pairs with an opening angle larger than 5 mrad were selected 
for the meson analysis. 

3 Neutral meson reconstruction 

Neutral pions and v\ mesons are reconstructed as excess yields, visible as peaks at their respective rest 
mass, above the combinatorial background in the two-photon invariant mass spectrum. Invariant mass 
spectra demonstrating the 7T° and r\ mesons peak in some selected p t slices are shown in Fig [I] by the 
histogram. The background is determined by mixing photon pairs from different events and is normalized 




Fig. 1: Invariant mass spectra in selected p t slices in PCM (left) and PHOS (right) in the n and r/ meson mass 
regions. The histogram and the bullets show the data before and after background subtraction, respectively. The 
curve is a fit to the invariant mass spectrum after background subttaction. 

to the same event background at the right side of the meson peaks. A residual correlated background 
is further subtracted using a linear or second order polynomial fit. The invariant mass spectrum after 
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background subtraction, depicted by bullets in Figjlj was fitted to obtain the n and 77 peak parameters 
(a curve). The number of reconstructed % s (t?s) is obtained in each p t bin by integrating the background 
subtracted peak within 3 standard deviations around the mean value of the 71° (T]) peak position in the 
case of PHOS. In the PCM measurement the integration windows were chosen to be asymmetric (/%o- 
0.035 GeV/c 2 , m^+0.010 GeV/c 2 ) and 0^-0.047 GeV/c 2 , ra,, +0.023 GeV/c 2 ) to take into account the 
left side tail of the meson peaks due to bremsstrahlung. For the same reason in the case of PCM the 
full width at half maximum (FWHM) instead of the Gaussian width of the peak was used. We vary 
the normalization and integration windows to estimate the related systematic uncertainties. The peak 
position and width from the two analyses compared to Monte Carlo simulations are shown in Fig. [2] as a 
function of p t . 




PCM • FWHM/2.36 
0— a) PHOS ■ □ a — 




1 10 p ( (GeV/c) 



Fig. 2: Reconstructed n° peak width (a) and position (b) as a function of p t in pp collisions at y/s = 7 TeV in 
PHOS and in the photon conversion method (PCM) compared to Monte Carlo simulations. The horizontal line in 
(b) indicates the nominal K° mass. 



The reconstruction efficiency £ and acceptance A are calculated in Monte Carlo simulations tuned to 
reproduce the detector response. In the PHOS case, the tuning included a 4.5% energy non-linearity 
observed in real data at E < 1 GeV and not reproduced by the GEANT simulations and an additional 6% 
channel-by-channel decalibration. In the PCM case, an additional smearing in each momentum compo- 



nent given by a = y o$ + of ■ p 2 with Co = 0.011 GeV/c and 0\ = 0.007 was necessary to reproduce 

the measured width of the 71° peak. PYTHIA [ 10] and PHOJET ifTTTl event generators and single particle 
simulations were used as input. The small photon conversion probability of about 8.5%, compensated by 
the large TPC acceptance, translates into £ - A of about 2 x 10~ 3 at p t > 1 GeV/c and decreases at lower 
p t due to the decrease of the efficiency of soft electron reconstruction and conversion probability. In the 
PHOS case, the acceptance A is zero for p t < 0.4 GeV/c, £ -A increases with p t and saturates at about 
2.0 x 10~ 2 at p t > 15 GeV/c. At high p t > 25 GeV/c the efficiency decreases due to cluster merging. 

The invariant differential cross section of 71° and 77 meson production were calculated as 

dp 3 2% Advents PtSABr Ay Ap t ' 

where cJmb r * s tne interaction cross section for the MBor trigger for pp collisions at \fs = 0.9 TeV or 
y/s = 7 TeV, Advents is the number of MBor events, p t is the transverse momentum within the bin to 
which the cross section has been assigned after the correction for the finite bin width Ap t (see below), 
Br is the branching ratio of the k° (17) meson to the two 7 decay channel and Af^ ™ is the number 
of reconstructed n° (17) mesons in a given Ay and Ap t bin. Finally, the invariant cross sections were 
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corrected for the finite p t bin width following the prescription in [24], keeping the y values equal to 
the bin averages and calculating the p x position at which the differential cross section coincides with 
the bin average. The Tsallis fit (see below) was used for the correction. Secondary 7T°'s from weak 
decays or hadronic interactions in the detector material are subtracted using Monte Carlo simulations. 
The contribution from K° decays is scaled using the measured K° spectrum at yfs = 0.9 TeV [30] or the 
charged kaon spectra at ^fs = 7 TeV |[3T| . The measured 71° and tj meson spectra at the center-of-mass 
energy of = 7 TeV cover a p t range from 0.3 to 25 GeV/c and from 0.4 to 15 GeV/c, respectively; the 
n° spectra at yfs = 0.9 TeV cover a p t range from 0.4 to 7 GeV/c. 





PHOS 




A = 1-1 GeV/c 


pt = 7.5 GeV/c 


Yield extr. 


±2.1 


±2.5 


Non-linearity 


±9.0 


±1.5 


Conversion 


±3.5 


±3.5 


Absolute energy scale 


±0.7 


±1.0 


Acceptance 


±1.0 


±1.0 


Calibration and alignment 


±7.0 


±3.0 


Pileup 


±0.8 


±0.8 


Total 


±12.5% 


±6.0% 




PCM 




p t =l.l GeV/c 


p t =7.5 GeV/c 


Material Budget 


±9.0 


±9.0 


Yield extraction 


±0.6 


±4.9 


PID 


±0.1 


±5.4 


x 2 (r) 


±0.3 


±6.2 


Reconstruction e 


±1.9 


±4.9 


Total 


±9.2% 


±14.0% 



Table 2: Summary of the relative systematic errors for the PHOS and the PCM analyses. 



A summary of the systematic uncertainties is shown in Table [2] for two different p t values. In PHOS, 
the significant source of systematic errors both at low and high p t is the raw yield extraction. It was 
estimated by varying the fitting range and the assumption about the shape of the background around the 
peak. The uncertainty related to the non-linearity of PHOS which dominates at low p t was estimated 
by introducing different non-linearities into the MC simulations under the condition that the simulated 
/^-dependence of the 71° peak position and peak width is still consistent with data. The uncertainties 
on the calibration and alignment were estimated in Monte Carlo simulations by varying the calibration 
parameters and the relative module positions within the expected tolerances. The uncertainty related to 
the pileup event rejection was evaluated in data by estimating the fraction of unidentified pileup events 
by extrapolating the distance and the number of contributing tracks of found pileup vertices to zero. 
The uncertainty of the conversion probability was estimated comparing measurements without magnetic 
field to the standard measurements with magnetic field. In the measurements with converted photons, 
the main sources of systematic errors are the knowledge of the material budget (dominant at low p t ), 
raw yield extraction, PID, the photon % 2 cut and reconstruction efficiency. The contribution from the 
raw yield extraction was estimated by changing the normalization range, the integration window, and 
the combinatorial background evaluation. The PID, photon % 2 {y) cut and reconstruction efficiency was 
estimated by evaluating stability of the results after changing the cut values. 
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Meson 


TeV 


A 


C 

(MeV/c 2 ) 


n 


7T° 


0.9 


1.5 + 0.3 


132±15 


7.8±0.5 


7T° 


7 


2.40±0.15 


139±4 


6.88±0.07 


n 


7 


0.21±0.03 


229 ±21 


7.0±0.5 



Table 3: Fit parameters of the Tsallis parametrisation (|2|) to the combined invariant production yields of n and 
rj mesons for inelastic events. The uncertainty on the parameter A due to the spectra normalization of and 
^35% at yfs = 900 GeV, and 7 TeV respectively, is not included. 



4 Results and comparison with pQCD 

The combined spectrum is calculated as a weighted average using statistical and systematic errors of the 
individual analyses |[25l . The combined production cross sections are shown in Fig.[3]a). The combined 
spectra including statistical and systematic errors are fitted with the Tsallis function [26 ] 

£ d%_ a^ EL i c-(n-l)(n-2) / | wit- m \~ " 
dp 3 2% nC [nC + m(n — 2)] y nC J ' 

where the fit parameters are A, C and n, a pp is the proton-proton inelastic cross section, m is the meson 
rest mass and m t = \/m 2 + pl is the transverse mass. The fit parameters are shown in Table [3] The 
property of the Tsallis function Q is such that the parameter A is equal to the integral of this function 
over p t from to infinity, A = dN/dy, and thus can be used as an estimation of the total yield at y = 
per inelastic pp collision. The additional uncertainty on the parameter A due to the spectra normalization 
of _{\% and 13 5'% at yfs = 900 GeV, and 7 TeV respectively, is not included. The found parameters of 
the Tsallis function for 71° production spectrum in pp collisions at y/s = 900 GeV are in agreement with 
those for the n + + %~ spectra measured by the ALICE collaboration at the same energy ll27ll . 

The ratio of the data points of the two methods to the combined fit, shown in Fig. |4j illustrates the 
consistency between the two measurements. 

We compare our results with Next-to-Leading Order (NLO) pQCD calculations using the PDF CTEQ6M5 
and DSS 7T° H, BKK 7T° JH and AESSS 17 [91 NLO fragmentation functions, see Fig.|]a). The data 
and NLO predictions are compared via a ratio with the fit to the measured cross section. This is shown in 
the bottom panels (b), (c) and (d) in Fig. [3] In the NLO calculations the factorization, renormalization and 
fragmentation scales are chosen to have the same value given by \i. The uncertainty in the inelastic pp 
cross section is represented by the full boxes at unity. At \f~s = 0.9 TeV the NLO calculations at pL = 1 p t 
describe the measured 71° data well, while at = 7 TeV the higher scale {\l=2 p t ) and a different set of 
fragmentation functions are required for a description of the data. However, the latter parameter set does 
not provide a good description of the low energy data. In any case, the NLO pQCD calculations show 
a harder slope compared to the measured results. Using the INCNLO program |[28l . we tested differ- 
ent parton distribution functions (CTEQ5M, CTEQ6M, MRS99) and different fragmentation functions 
(BKK, KKP, DSS) and found a similar result: pQCD predicts harder slopes, and variation of PDFs and 
FFs does not change the shape, but results mainly in the variation of the absolute cross section. A similar 
trend is observed for the v\ meson (a higher scale /I = 2p t is required), although the discrepancy is less 
significant due to the larger error bars and smaller p t reach. 

The ratio 77 /tt° is shown in Fig. [5] It has the advantage that systematic uncertainties in the measurement 
partially cancel. This is also the case for the NLO pQCD calculation, where in particular the influence of 
the PDF is reduced in the ratio. Here, predictions that failed to reproduce the measured 71° and v\ cross 
section are able to reproduce the 17/71° ratio. 
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Fig. 3: a) Differential invariant cross section of n production in pp collisions at 



7 TeV (circles) and 0.9 TeV 



(squares) and of rj meson production at = 7 TeV (stars). The lines and the boxes represent the statistical 
and systematic error of the combined measurement respectively. The uncertainty on the pp cross section is not 
included. NLO pQCD calculations using the CTEQ6M5 PDF and the DSS (AESS for rj mesons) FF for three 
scales jj. — 0.5p t , lp t and 2p t are shown. Dotted lines in panels b) and c) correspond to the ratios using the BKK 
FF. Ratio of the NLO calculations to the data parametrisations are shown in panels b), c) and d). The full boxes 
represent the uncertainty on the pp cross sections. 
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Fig. 4: Ratio of the two independent n meson measurements to the fit of the combined normalized invariant 
production cross section of n° mesons in pp collisions at \/s = 7 TeV. 
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Fig. 5: 77/ 7r° ratio measured in pp collisions at y/s = 7 TeV compared to NLO pQCD predictions. 



5 Conclusion 

In summary, the invariant differential cross sections for inclusive tt° production in pp collisions at y/s = 
7 TeV and 0.9 TeV and for tj meson production at 7 TeV have been measured in a wide p t range taking 
advantage of two independent methods available in the ALICE experiment at the LHC. NLO pQCD 
calculations cannot provide a consistent description of measured data at both beam energies. State-of- 
the-art calculations describe the data at 0.9 TeV and 0.2 TeV ll32l . however this is not the case at 7 TeV, 
where the calculations overestimate the cross sections and exhibit a different slope compared to the data. 
Thus, this measurement provides an important input for the tuning of pQCD calculations and represents 
crucial reference data for the measurement of the nuclear modification factor Raa of the 71° production 
in heavy-ion collisions at the LHC. Furthermore, the NLO predictions for the tj mesons using the newest 
fragmentation functions require a value jj. = 2p t in order to get closer to the experimental results. 
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